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ABSTRACT 

Context. After several generations of interferometers in radioastronomy, wide-field imaging at high angular resolution 
is today a major goal for trying to match optical wide-field performances. 

Aims. All the radio-interferometric, wide-field imaging methods currently belong to the mosaicking family. Based on a 
^ . 30 years old, original idea from Ekers & Rots, we aim at proposing an alternate formalism. 

Methods. Starting from their ideal case, we successively evaluate the impact of the standard ingredients of interferometric 
■^L ' imaging, i. e. the sampling function, the visibility gridding, the data weighting, and the processing of the short spacings 

either from single-dish antennas or from heterogeneous arrays. After a comparison with standard nonlinear mosaicking, 
we assess the compatibility of the proposed processing with 1) a method of dealing with the effect of celestial projection 
and 2) the elongation of the primary beam along the scanning direction when using the on-the-fly observing mode. 
Results. The dirty image resulting from the proposed scheme can be expressed as a convolution of the sky brightness 
distribution with a set of wide-field dirty beams varying with the sky coordinates. The wide-field dirty beams are locally 
shift-invariant as they do not depend strongly on position on the sky; Their shapes vary on angular scales typically 
larger or equal to the primary beamwidth. A comparison with standard nonlinear mosaicking shows that both processing 
(— I ' schemes are not mathematically equivalent, though they both recover the sky brightness. In particular, the weighting 

scheme is very different in both methods. Moreover, the proposed scheme naturally processes the short spacings from 
I . both single-dish antennas and heterogeneous arrays. Finally, the sky gridding of the measured visibilities, required by 

O ' the proposed scheme, may potentially save large amounts of hard-disk space and cpu processing power over mosaicking 

^ ' when handling data sets acquired with the on-the-fly observing mode. 

C/3 , Conclusions. We propose to call this promising family of imaging methods wide-field synthesis because it explicitly 

■ synthesizes visibilities at a much finer spatial frequency resolution than the one set by the diameter of the interferometer 

antennas. 

Key words, techniques: wide-field - techniques: interferometric - methods: data analysis - methods: image processing 
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[ 1. Introduction The ideal measurement equation of interferometric 

1 . . . . wide-field imaging is 

, The instantaneous field of view of an interferometer is nat- 

' urally limited by the primary beam size of the individual p 

O ' antennas. For the ALMA 12m-antennas, this field of view V{up,as) ^ / i3(ap - as) /(ap) e"*^'""""" da^, (1) 

O is - 9" at 690 GHz and - 27" at 230 GHz. The astro- -'"p 

' physical sources in the (sub)-millimeter domain are often 

j> ! much larger than this, but still structured on much smaller where V is the visibility function of 1) Up (the spatial 

K> ' ^^S^^'"^^ scales. Interferometric wide-field techniques enable frequency with respect to the fixed phase center) and 2) 

^ ; us to fully image these sources at high angular resolution, (the scanned sky angle), / is the sky brightness, and 

^ : These techniques first require an observing mode that in B the antenna power patte rn or primary beam o f an an- 

■ - - ' one way or another scans the sky on spatial scales larger tenna of the interferometer (jThompson et al.lll98a chapter 

than the primary beam. The most common observing mode 2). For simplicity, 1) we assume that the primary beam 

in use today, known as stop-and-go mosaicking, consists in is independent of azimuth and elevation, and 2) we use 

repeatedly observing sky positions typically separated by one-dimensional notation without los s of generality. We do 

half the primary beam size. The improvement of the track- not deal with polarime try (see e.g. iHamaker et al.l 119961: 

ing behavior of modern antennas now leads astronomers to ISault et aLlll996al Il999f) because it is adds another level of 

consider on-the-fly observations, with the antennas slewing complexity over our first goal here, i.e. wide-field consid- 

continuously across the sky. The improvements in correla- erations. Several aspects make Eq. [T] peculiar with respect 

tor and receiver technologies are also leading to techniques to the ideal measurement equation for single-field observa- 

that could potentially sample the antenna focal planes with tions. First, the visibility is a function not only of the uv 

multi-beam receivers instead of the single-pixel receivers in- spatial frequency (up) but also of the scanned sky coordi- 

stalled on current interferometers. nate (as). Second, Eq. [T]is a mix between a Fourier trans- 

form and a convolution equation. It can be regarded, for 

Send offprint requests to: e-mail: pety@iram.fr example, as the Fourier transform along the ttp dimension 
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of the function, B{ap — as) I{ap), of the {ap,as) variables. 
But Eq. [T] can also be written as the convolution: 



V{up, as 



B{as - ap) I{ap, Up) dap, 



where B{as — ap) = B{ap — Ug) 
and I(ap,itp) EE /(ap)e-*27rapup^ 



(2) 

(3) 
(4) 



For each Up kept constant, V{up,as) is the convolution of 
B and T. Indeed, X{ap,Up = 0) = I{ap), so we derive 



V{up^{),as)= f B{c 



ap) I{ap) dap, 



(5) 



i.e. , the convolut i on equ ation for single-dish observations. 

ers fc Rots! (I1979D were the first to show that the 
measurement equation (Eq. [1]) enables us to recover spa- 
tial frequencies of the sky brightness at a much finer uv 
resolution than the uv resolution set by the diameter of 
th e interferometer ante nnas. Interestingly enough, the goal 
of lEkers fc Rots! ()1979[ ) was "just" to find a way to pro- 
duce the missir ig short sp acings of a multiplying interferom- 
eter. However, ICornwelll (|T988i ) realized that Ekers & Rots' 
scheme has a much stronger impact, because it exp lains why 
an interferometer is able to do wide-field imaging. ICornweUl 
(|198 8) also demonstrated that on-the-fiy scanning is not 
absolutely necessary to interferometric wide-field imaging. 
Indeed, the large-scale information can be retrieved in mo- 
saics of single-field observations, provided that the sampling 
of the single fields follows the sky-plane Nyquist sampling 
theorem. 

As a result, all the information about the sky brightness 
is coded in the visibility function. From a data-processing 
viewpoint, all the current radio- interferometr i c wide- field 
imaging methods (see, e.g . , Guct h et al.l 119 95': Sault et al.' 
1996bHCornwell e t all IT993; Bhatnagar & Cornwell ,2004; 
Bhatnagar et al.r 2008: Cotton & Uson 2008) belong to the 
mosaicking familjQ pioneered by iCornwelli (| 19881 ) . In this 



family, the processing starts with Fourier transforming 
V{up,as) along the Up dimension (i.e. at constant as) to 
produce a set of single-field dirty images before linearly 
combining them and forming a wide-field dirty image. In 
this paper, we propose an alternate processing, which starts 
with a Fourier transform of V{up,as) along the as dimen- 
sion (i.e. at constant Up). We show how this explicitely 
synthesizes the spatial frequencies needed to do wide-field 
imaging, which are linearly combined to form a "wide-field 
uv plane", i.e., one uw-plane containing all the spatial fre- 
quency information measured during the wide-field obser- 
vation. Inverse Fourier transform will produce a dirty im- 
age, which can then be deconvolved using standard meth- 
ods. The existence of two different ways to extract the wide- 
field information from the visibility function raises several 
questions: Are they equivalent? What are their relative 
merits? 

We thus aim at revisiting the mathematical foundations 
of wide-field imaging and deconvolution. Sections[2]to[7]pro- 
pose the new algorithm, which we call wide-field synthesis: 
Section [5] first defines the notations and it then lays out the 



^ In the rest of this paper, stop-and-go mosaicking refer to 
the observing mode, while mosaicking alone refer to the imaging 
family. 



basic concepts used throughout the paper. Section [3] states 
the steps needed to go beyond the Ekers & Rots scheme and 
explores the consequences of incomplete sampling of both 
the uv and sky planes. Section H] discusses the effects of 
gridding by convolution and regular resampling. Section [5] 
describes how to influence the dirty beam shapes and thus 
the deconvolution. Section [5] states how to introduce short 
spacings measured either from a single-dish antenna or from 
heterogeneous interferometers. Section [7] compares the pro- 
posed wide-field synthesis algorithm with standard non- 
linear mosaicking. Some detailed demonstrations are fac- 
tored out in Appendix |X] to enable an easier reading of 
the main paper, while ensuring that interested readers can 
follow the demonstrations. Appendices [B] and [O then ex- 
plain how the wide-field synthesis algorithm can cope with 
non-ideal effects: Appendix [B] discusses how at least one 
standard way to cope with sky projection problems is com- 
patible with the wide-field synthesis algorithm. Appendix [Cl 
explores the consequences of using the on-the-fly observ- 
ing mode. Finally, we assume good familiarity with single- 
field imaging in various places. We refer the reader to well- 
kno wn references: e.g. chapte r 6 of iThompson et aD (|l986f) 
and lSramek fc Schwab! (|1989t) . 

2. Notations and basic concepts 

2.1. Notations 

In this paper, we use the H racewe l (pool 's notation to 
display the relationship between a function I{a) and its 
direct Fourier transform I{u), i.e., 



lia) D I{u), 



(6) 



where {a,u) is the couple of Fourier conjugate variables. We 
also use the following sign conventions for the direct and 
inverse Fourier transforms 



I{u)= / /(a) e-*2^"" da 
and 



I{a) 



I{u)e 



(7) 



(8) 



As y is a function of two independent quantities {up and 
Us), the Fourier transform may be applied independently 
on each dimension, while the other dimension stays con- 
stant. Several additional conventions are used to express 
this. First, we introduce a specific notation to state that 
either the first or the second dimension stays constant: 

Kp(as) = V{up = constant, as), (9) 

and 

V"^{up) = V{up,as = constant). (10) 

Second, we use a bottom/top line to derive the notation 
of the Fourier transform along the first/second dimension 
from the notation of the original function. Third, on the 
Fourier transform sign [i.e. d), we explicitly state the di- 
mension along which the Fourier transform is computed. 
For instance, if Z? is a function of {ap,as), then the Fourier 
transform of D along the first dimension is expressed as 



D{ap,as) p D{up,as), 



(11) 



2 



J. Pety and N. Rodriguez-Fernandez: Revisiting the theory of interferometric wide-field synthesis 



Table 1. Definition of the symbols used to expose the wide- 
field synthesis formalism. 



Symbol & Definition Plane(s)" 



Ota 


Scanned angle 


sky 


Us 


Scanned spatial frequency 


uv 


Ctp 


Phased angle 


sky 


Up 


Phased spatial frequency 


uv 


I 


Sky brightness 


sky 


B 


Primary beam 


sky 


V 


Visibility function 


uv & sky 


s 


Sampling function 


uv & sky 


A 


Set of single-field dirty beams 


sky & sky 


D 


Set of wide-field dirty beams 


sky & sky 




Sky-plane weighting function 


sky & sky 


W 


un-plane weighting function {Q. D W) 


uv & uv 


Q 


Gridding function (= g 7) 


uv & sky 


9 


ui)-plane gridding function 


uv 


7 


Sky-plane gridding function 


sky 


-^dirty 


Wide-field dirty image 


sky 



Notes. Planes of definition of the associated symbols. 



while the Fourier transform of D along the second dimen- 
sion is expressed as 

as 

D{ap,as) D D{ap,Us). (12) 

We also use a more compact notation when doing the 
Fourier transform on both dimensions simultaneously, i.e., 

D{ap,as) ^^^^^^ D{up,u,). (13) 

Finally, the convolution of two functions G and V is noted 
and defined as 

{G*V}iu)= f Giu-v)V{v)dv. (14) 

J V 

For reference, Table [1] summarizes the definitions of the 
symbols used most throughout the paper. With the one- 
dimensional notation used throughout the paper, the num- 
ber of planes quoted directly gives the number of associated 
dimensions of the symbols. Generalization to images would 
require a doubling of the number of planes/dimensions. 
Table [5] defines the uv and angular scales that are relevant 
to wide-field interferometric imaging, and Fig. [1] sketches 
the different angular scales. Each angular scale (9) is re- 
lated to a uv scale (d) through 6 — 1/d, where 6 and d 
are measured in radians and in units of A (the wavelength 
of the observation). In the rest of the paper, we explicitely 
distinguish between ^prim = 1/rfprim, the angular scale as- 
sociated to the diameter of the interferometer antennas, 
and 0fwhmj the full width at half maximum of the primary 
beam. The relation between 6'prim and 6'twhm depends on 
the illumination of the receiver feed by the antenna optics. 
In r adio astronomy, w e typically have ^fwhm ~ 1.2 0prim (see 
e.g. iGoldsmitbl I1998L chapter 6). Finally, the notion of 
anti-aliasing scale (^aiias) is introduced and discussed in 
Sect. 1121 

2.2. Basic concepts 

Figure [5] illustrates the principles underlying 1) the setup 
to get interferometric wide-field observations and 2) our 



Table 2. Definition of the uv and sky scales relevant to 
wide-field interferometric imaging. 



Symbol Definition 

[A,rad]'' Conjugate uv and angular scale 

dmax,6'syn Maxlmum baseline length & Synthesized beam 

dprim, Sprim Autenua diameter & Primary beamwidth 

daiiaa , ^aiias Minimum image size for tolerable aliasing 

rffloid , ^floid Targeted field of view 

rfimage, ^image Final Image size 



Notes. The chosen units (radians for 9 and wavelength for 
d) imply that the conjugate scales are linked through 9 = 1/d, 
instead of the usual 9 = X/d. 







Synthesized beamwidth 

(e.^ = i/d„„) 






^^ I'rimary beamwidth 

(^fwhm = l-'^^prim = ^'^/'^prim) 






/ \ Minimum dirty beam image size 

(e-g- 9aiias = 2.2 Sf^iiJ for 
V J tolerable abasing (< 10 ) 






Targeted field of view (0fieid) 




Image size {9^^^^^ = efidd+Saiias) 



Fig. 1. Visualization of the different angular scales rele- 
vant to wide-field interferometric imaging. The notion of 
anti-aliasing scale (^aiias) is introduced and discussed in 
Sect.lU 

proposition to process them. For simplicity, we display the 
minimum possible complexity without loss of generality. 
The top row displays the sky plane. The middle row rep- 
resents the 4-dimensional measurement space at different 
stages of the processing. As it is difficult to display a 4- 
dimensional space on a sheet of paper, the bottom row 
shows 2-dimensional cuts of the measurement space at the 
same processing stages. 

2.2.1. Observation setup and measurement space 

Panel a) displays the sky region for which we aim for es- 
timating the sky brigthness, I (a). The field of view of an 
interferometer observing in a given direction of the sky has 
a typical size set by the primary beam shape. In our ex- 
ample, this is illustrated by any of the circles whose diam- 
eter is 0prim- As we aim at observing a wider field of view, 
e.g. Afield, the interferometer needs to scan the targeted sky 
field. We assume that we scan through stop-and-go mo- 
saicking, ending up with a 7-field mosaic. 

After calibration, the output of the interferometer is 
a visibility function, V{up,as), whose relation to the sky 
brightness is given by the measurement equation (Eq. [T|). 
Panel b.l) shows the measurement space as a mosaic of 
single-field uv planes: The uv plane coverage of each single- 
field observation is displayed as a blue sub-panel at the 
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e) 

efield/2 




-efield/2 




I(a,|S) 



Scanned 
interferometric 
observation 



ide — field synthesis 



Iwide(«.^) 



1. Fourier 
transform 

2. Deconvolution 



Fourier transform 1. Shift 
along and /Sg , 2. Average 
V(Up,Vp,Oi„^J > V(Up,Vp,U3,V3) > <V>(u,v) = Idirty(u,v) 




1 "s2 "s3 




d. 





1 






















f 
> 








II 

> 










1 







Up [A] 



'J = Up+U3 [A] 



b.2 



C.2 



li.ld/2 -"pr,m/2 



, [rod] 

(ield/2 













\ 






\ 








N 




=u 






1 


\\ 






\ 


-dfieia 














\ \ 






\ 





d.2 



"="p+". w 

I I > 



Fig. 2. Illustration of the principles of wide-field synthesis, which enables us to image wide-field interferometric observa- 
tions. The top row displays the sky plane. The middle row displays the 4-dimensional visibility space and the bottom row 
shows 2-dimensional cuts of this space at different stages of the processing. In panels b) to d), the scanned dimensions 
{as and Ug) are displayed in blue while the phased spatial scale dimensions (up) are displayed in red and the spatial 
scale dimensions (it) of the final wide-field uv plane are displayed in black. The grey zones of panels b.2) and c.2) show 
the regions of the visibility space without measurements (missing short-spacings). In detail, panel a) shows a possible 
scanning strategy of the sky to measure the unknown brightness distribution at high angular resolution: For simplicity 
it is here just a 7-field mosaic. Panel b.l) and b.2) sketch the space of measured visibilities: The uv plane at each of 
the 7 measured sky positions is displayed as a blue square box in panel b.l) and a blue vertical line in panel b.2). For 
simplicity, only 6 visibilities are plotted in panel b.l). Panels c.l) and c.2) sketch the space of synthesized visibilities after 
Fourier transform of the measured visibilities along the scanned coordinate (as): At each measured spatial frequency Up 
(displayed on the blue axes) is associated one space of synthesized wide-field spatial frequencies displayed as one of the 
red squares in panel c.l) and the red vertical lines in panel c.2). The wide-field spatial scales are synthesized 1) on a 
grid whose cell size is related to the total field of view of the observation and 2) only inside circles whose radius is the 
primary diameter of the interferometer antennas. Panels d.l) and d.2) display the final, wide-field uv plane. This plane is 
built by application of the shift-and-average operator along the black lines on panel c.2), lines that display the region of 
constant u spatial frequency in the {up,Us) space. Standard inverse Fourier transform and deconvolution methods then 
produce a wide-field distribution of sky brightnesses as shown in panel e). 
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sky position where it has been measured and which is fea- 
tured by the red axes. We assume 1) that the interfer- 
ometer has only 3 antennas and 2) that only a single in- 
tegration is observed per sky position. This implies only 
6 visibilities per single-field uv plane. In panel b.2), the 
uv planes at constant are displayed as the blue verti- 
cal lines. The measured spatial frequencies belong to the 
[-dmax, -rfmin] and [+dmin , +'imax] ranges, where dmin and 
c^max are respectively the shortest and longest measured 
baseline length, dmin is related to the minimum tolerable 
distance between two antennas to avoid collision. Here, we 
chose dmin ~ l-5dpriin- The grey zone between — dmin, and 
-|-ciniin displays the missing short spacings. 

2.2.2. Processing by explicit synthesis of the wide-field 
spatial frequencies 

All the information about the sky brightness, I (a), is some- 
how coded in the visibility function, V{up,as). The high 
spatial frequencies (from dmin to dmax) are clearly coded 
along the Up dimension. The uncertainty relation between 
Fourier conjugate quantities also implies that the typical 
spatial frequency resolution along the Up dimension is only 
dprim because the field of view of a single pointing has a 
typical size of 6'priin- However, wide-field imaging implies 
measuring all the spatial frequencies with a finer resolu- 
tion, dfieid = 1/^fieid- The missiug information must then 
be hidden in the as dimension. 

In Sect. [31 we show that Fourier transforming the mea- 
sured visibilities along the dimension (i.e. at constant 
Up) can synthesize the missing spatial frequencies, because 
the Us dimension is sampled from —Oucm/'^ to -f6'ficid/2, 
implying a typical spatial-frequency resolution of the Us 
dimension equal to dfidd- Conversely, the dimension is 
probed by the primary beams with a typical angular res- 
olution of 0piim, implying that the Us spatial frequencies 
will only be synthesized inside the [— dprim, +dprim] range. 
Panels c.l) and c.2) illustrate the effects of the Fourier 
transform of V{up,Us) along the as dimension, in 4 and 2 
dimensions, respectively. The red subpanels or vertical lines 
display the Us spatial frequencies around each constant Up 
spatial frequency. 

In panels d.l) and d.2) {i.e. after the Fourier trans- 
form along the as dimension), V{up,Us) contains all the 
measured information about the sky brightness in a spatial 
frequency space. However, the information is ordered in a 
strange and redundant way. Indeed, we show that V{up,Us) 
is linearly related to I{up+Us). To first order, the informa- 
tion about a given spatial frequency u is stored in all the 
values of V{up^Us) which verifies u ^ Up + Ug (black lines 
on panel c.2). 

A shift operation will reorder the spatial scale infor- 
mation and averaging will compress the redundancy (illus- 
trated by the halving of the number of the space dimen- 
sions). The use of a shift- and- average operator thus pro- 
duces a final uv plane containing all the spatial scale in- 
formation to image a wide field in an intuitive form. We 
thus call this space the wide-field uv plane. Panels d.l) and 
d.2) display this space, where the minimum relevant spa- 
tial frequency is related to the total field of view, while the 
maximum one is related to the interferometer resolution. 

Sections [3] and ^ show that applying the shift-and- 
average operator to V produces the Fourier transform of a 
dirty image, which is a local convolution of the sky bright- 



ness by a slowly varying dirty beam. As a result, inverse 
Fourier transform of (V') and deconvolution methods will 
produce a wide-field distribution of sky brightness as shown 
in panel e) at the top right of Fig. [51 

3. Beyond the Ekers & Rots scheme 

In the real world, the visibility function is not only sam- 
pled, but this sampling is incomplete for two main reasons. 
1) The instrument has a finite spatial resolution, and the 
scanning of the sky is limited, implying that the sampling 
in both planes has a finite support. 2) The uv coverage 
and the sky-scanning coverage can have holes caused either 
by intrinsic limitations (e.g. lack of short spacings or small 
number of baselines) or by acquisition problems (implying 
data flagging). The incomplete sampling makes the mathe- 
matics on the general case complex. We thus start with the 
ideal case where we assume that the visibility function is 
continuously sampled along the Up and as dimension. We 
then look at the general case. 

3.1. Ideal case: Infinite, continuous sampling 

Starti ng from the measurement equation [1] lEkers &: Rots! 
()1979D first demonstrated (see Sect. [XT]) thaci 

y{Up,Us), VupiUs) = B{~Us)7{Up+Us). (15) 

For each constant Up spatial frequency, the Fourier trans- 
form thus synthesizes a function, Vup{us), which is sim- 
ply related to I{up+Us), the Fourier components of the 
sky brightness around Up. V{up,Us) is only defined in the 
[— dprim, +dprim] interval along the Us dimension because 
B{-Us) is itself only defined inside this interval, since B{- 
Us) is the autocorrelation of the antenna illumination. 

We search to derive a single estimate of the Fourier com- 
ponents I{u) of the sky brightness. Equation [TSl indicates 
that the fraction V{up,Us)/ B{—Us) gives us an estimate of 
I{u) for each couple {up,Us) that satisfies u — Up + Us- 
However, the information about / is strangely ordered. 
There are two possible ways to look at this ordering. 1) 
Starting from the measurement space, the Ekers & Rots 
scheme synthesizes frequencies around each Up measure in- 
side the interval [up — dprim, Up + dprim] at the dfioid spatial 
frequency resolution. 2) Starting from our goal, we want to 
estimate / at a given spatial frequency u with a dfieid spa- 
tial frequency resolution. We thus search for all the couples 
{up, Us) satisfying u = Up + Us, which are displayed in panel 
c.2) of Fig.[21as the diagonal black lines. It immediately re- 
sults that 1) there are several estimates of / for each spatial 
frequency u and 2) the number of estimates varies with u. 
We can average them to get a better estimate of I{u). 

^ The convolution theorem, which states that the Fourier 
transform of the convolution of two functions is the product of 
the Fourier transform of both individual functions, is a special 
case for Eq. 1151 It can be recovered by setting Up = 0. Indeed, 
as already mentioned in the introduction, the ideal measure- 
ment equation [1] can be interpreted as a convolution with an 
additional phase term. By Fourier transforming along the Qa 
dimension, the convolution translates into a product of Fourier 
transforms B and /, while the phase term translates into a shift 
of coordinates: Up + Us. 
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This last viewpoint thus suggests averaging in the 3.1.2. Summary and interpretation 



{up,Us) space along linepaths defined hy u = Up 



Such 



an operator can mathematically be defined as 
{F){u)=ll 5[u - (up + Us)] W{up,Us) F{up, Us) dupdus, (16) 



where F is the function to be averaged and is a normal- 
ized weighting function. Using the properties of the Dirac 
function, we can reduce the double integral to 



(FKu) 



W{up, u — Up) F{up, u — Up) dup 



(17) 



In this equation, we easily recognize a shift- and- average op- 
erator. The normalized weighting function plays a critical 
role in the following formalism, and we propose clever ways 
to define W in Sect. [5] In the ideal case studied here, W 
can be defined as 



W{up,Us) = l/(2%/2dpri„i) for 
W{up,Us) = for other values of ' 



^prim 



primj ; 



In other words, we have just normalized the integral by the 
constant length (2-\/2 dprim) of the averaging linepath. 

3.1.1. Wide-field dirty image, dirty beam and image-plane 
measurement equation 



Section 13.21 shows that the incomplete sky and uv sam- 
pling_forbid us_to apply the shift-and-average operator to 
the V[up, Us)/ B{—Us) function. To guide us in this general 
case, we thus explore the consequences of applying this op- 
erator to V in the ideal case. It is easy to demonstrate that 
the result is the Fourier transform of a dirty image, i.e.. 



U,.ty{u)^{V){u). 



(18) 



Indeed, substituting (y){u) with the help of Eos. [TTlandfTSl 
and taking the inverse Fourier transform, we get 



' dirty 



{u)=D{u)I{u), 



with D{u) 



W{up, u — Up) B{up — u) dup. 



(19) 



(20) 



Here, /dirty conforms to the usual idea of dirty image, i.e.. 
the convolution of a dirty beam by the sky brightness: 



/dirty (a) = {£>★/} (a) , 



(21) 



In contrast to the usual situation for single-field observa- 
tions, the mix between a Fourier transform and a convo- 
lution of Eq. [TJ associated with the specific processing 
changes the image-plane measurement equation from a con- 
volution of a dirty beam with the product B I to a convo- 
lution of a dirty beam with /. The dependency on the pri- 
mary beam is still there. It is just transferred from a product 
of the sky brightness distribution into the definition of the 
dirty beam. 



I.e. direct Fourier transform along the Os dimension, shift- 
and-average to define a final wide-field uv plane, and inverse 
Fourier transform. 



In summary, a theoretical implementation of wide-field syn- 
thesis implies 

1. the possibility of Fourier transforming the visibility 
function along the as dimension (i.e. at constant Up), 
which gives us a set of synthesized uv planes; 

2. the possibility of shifting-and-averaging these synthe- 
sized uv planes to build the final, wide-field uv plane 
containing all the available information. 

Using those tools, we are able to write the wide-field image- 
plane measurement equation as a convolution of a wide-field 
dirty beam [D) by the sky brightness (/), i.e.. 



J dirty 



(a) 



D{a- a')I{a')da' . 



(22) 



We can write a convolution equation in this ideal case be- 
cause the wide-field response of the instrument is shift- 
invariant; i.e., D only depends on differences of the sky 
coordinates. 

It is well-known that for a single-field observation, the 
dirty beam is the inverse Fourier transform of the sampling 
function. The shape of this sampling function is due to the 
combination of aperture synthesis (the interferometer an- 
tennas give a limited number of independent baselines) and 
Earth-rotation synthesis (the rotation of the Earth changes 
the projection of the physical baselines onto the plane per- 
pendicular to the instantaneous line of sight). By analyzing 
via a Fourier transform, the evolution of the visibility func- 
tion with the sky position, the Ekers & Rots scheme syn- 
thesizes visibilities at spatial frequencies needed to image a 
larger field of view than the interferometer primary beam. 
We thus propose to call this specific processing: wide-field 
synthesis. 

3.2. General case: Incomplete sampling 

Reality imposes limitations on the synthesis of spatial fre- 
quencies. Indeed, we have already stated that the visibility 
function is incompletely sampled both in the uv and sky 
planes. To take the sampling effects into account, we in- 
troduce the sampling function S{up,as), which is a sum of 
Dirac functions at measured positiontQ. The sampling func- 
tion cannot be factored into the product of two functions, 
each only acting on one plane. Indeed, the Earth rotation 
happening during the source scanning implies a coupling of 
both dimensions of the sampling function. In other words, 
the uv coverage will vary with the scanned sky coordinate. 
This leads us to a shift-dependent situation, precluding us 
from writing the wide-field image-plane measurement equa- 
tion as a true convolution. We nevertheless search for a 
wide-field image-plane measurement equation as close as 
possible to a convolution because all the inversion meth- 
ods devised in the past three decades in radioastronomy 
are tuned to deconvolve images. The simplest mathemati- 
cal way to generalize Eq. [21] to a shift-dependent situation 
is to write it as 



/dirty (a) 



D{a~ a',a)I{a')da'. 



(23) 



* Loosely speaking, the sampling function can be thought as 
a function whose value is 1 where there is a measure and 
elsewhere. 
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In this section, we show how the linear character of the 
imaging process allows us to do this. Section 13.2.11 derives 
the impact of incomplete sampling on the Ekers & Rots 
equation, and Sect. 13.2.21 derives the wide-field measure- 
ment equation in the uv plane. Section r3.2.3l interprets these 
results. 



3.2.1. Effect on the Ekers &l Rots equation 

The sampled visibility function, SV, is defined as the prod- 
uct of S and V and SV its Fourier transform along as, i.e., 



SV{up, as) = S{up, as) V{up, as), 
and SV{up,as) D SV{up,Us). 



(24) 
(25) 



Because SVu^, is the product of two functions of a^, we 
can use the convolution theorem to show that SVup is the 
convolution of Su„ by Vu„, i.e.. 



SVu^iUs) 



Su^{us ~ u^)Vup{u's) du's- 



(26) 



By replacing Vup with the help of the Ekers & Rots relation 
(Eq. [151) J derive 



SVupiUs) 



S^^ (us - <) i?(-<) liup + u's) d<. (27) 



As B is bounded inside the [— dprim, +c?prim] interval, 
SVupius) is a local average, weighted by Sup{us — 
u'j,) B{—u'^), of I{up+u'^) around the Up spatial frequency. 

As expected, we recover Eq. [TS] for the ideal case 
(i.e., infinite, continuous visibility function) because then 
Supius — u's) = d{us — u's). A more interesting case arises 
when the visibility function is continuously sampled over a 
limited sky field of view, i.e.. 



\/up, S{up,as)^l if la^l < 6'ficid/2, 
Vup, S{up,as)=0 if la^l > 6'ficid/2. 

After Fourier transform this gives 

1 . 



Vup, S{up,Us) 



■ smc 



icld 



Afield 



(28) 
(29) 



(30) 



In this case, the local average of the sky brightness Fourier 
components happens on a typical uv scale equal to cZfieid- 
However, the sine function is known to decay only slowly. 
Some observing strategy (e.g. quickly observing outside the 
edges of the targeted field of view to provide a bandguard) 
could be considered to apodize the sky-plane dependence of 
the sampling function, resulting in faster decaying S func- 
tions, hence in less mixing of the wide-field spatial frequen- 
cies. 



3.2.2. ww-plane wide-field measurement equation 

Because we aim at estimating the Fourier component of /, 
we introduce the following change of variables u' = Up + u's 
and du' — du'^, to derive 

SVuj,{us) = I Supiup -I- Ms - u') B(up - u')l{u') du' . (31) 



We then shift-and-average SV{up,Us) to build the Fourier 
transform of a wide-field dirty image 



dirty (m) = {SV){u), with U— Up+Us 



(32) 



Substituting the shift-and-average operator by its definition 
and using Eq. [3T]to replace SVu^ius), we derive 

^dirty(M) 



W{up,u — Up) S{up, u — u') B{up — u') I{u') dupdu' .{33) 



This uw-plane wide-field measurement equation can be 
written as 

/dirty (m) = [ mu', U - u') 7{u') du' , (34) 
J u' 

if we enforce the following equality 

D.{u' , u—u') = W{up, u~Up) S{up, u~u') B{up—u') dup.(35) 

J Up 

This is one way to define D_, which is convenient though 
unusual. It is implicit in this definition that we need to 
make a change of variable {u" — u ^ u') to derive 

Qiu', u") = I W{up, u'+u"~up) S{up, u") B{up-u') dup.{36) 

J Up 

In the following, we use either one or the other definition 
of D_, depending on convenience. 

3.2.3. Interpretation 



Appendix IA.2I demonstrates that the image and uw-plane 
wide-field measurement equations (Eqs. [221 and IMl) are 
equivalent if 



{ap,as) 

D{ap,as) ^^^^^^^ D{up,Us). 



(37) 



The image-plane wide-field measurement equation (Eq. (231) 
can be written as 

/dirty(a) */}(«). (38) 

Its interpretation is straightforward: The sky brightness dis- 
tribution is convolved with a dirty beam, D{a' ,a"), which 
varies with the sky coordinate a" . This raises the question 
of the rate of change of the dirty beam with the sky coor- 
dinate. This question is addressed in Sects. and [51 

4. Gridding by convolution and regular resampling 

We want to Fourier transform the raw visibilities along the 
sky dimension (as) at some constant value in the Up dimen- 
sion. The raw data, however, is sampled on an irregular grid 
in both the uv and sky planes. We need to grid the mea- 
sured visibilities in both the uv and the sky planes before 
Fourier transformation for different reasons. First, the grid- 
ding in the uv plane will handle the variation in the spatial 
frequency as the sky is scanned, i.e., the difficulty and per- 
haps the impossibility of Fourier-transforming at a com- 
pletely constant Up value. Second, the gridding along the 
sky dimension allows the use of Fast Fourier Transforms. 
As usual, we grid through convolution and regular resam- 
pling. 



7 



J. Pety and N. Rodriguez-Fernandez: Revisiting the theory of interferometric wide-field synthesis 



4.1. Convolution 
4.1.1. Definitions 

We first define a gridding kernel that depends on both 
dimensions, Q{u,as). This gridding kernel can be chosen 
as the product of two functions, simplifying the following 
demonstrations: 



(39) 



We then define the sampled visibility function gridded in 
both the uv and sky planes as 

SV^iup,as) = {g*SV}iup,as) (40) 

g{up - Up) j{a, - a'J SV{u'p, a'^) du'^da'^. (41) 

Finally, when assessing the impact of the gridding on the 
measurement equation 1341 a new function, 



E(up,as,a") = S{up,as) B{a'!. - a^), 



(42) 



and its Fourier transforms naturally appear in the equa- 
tions. Defining the following Fourier transform relationships 



E(up,as,a") J?, S(up,Ms,a"), 
and 

E(up,M^,a") D E(up,U5,-u"), 
we easily derive 

E(up,a„0 = 5(7/p,a,)B«)e-^2.<a, 
and 



(43) 
(44) 

(45) 
(46) 



Using these notations, we have before gridding, 

SViup, as) = [ E(up, as,ap) I{ap) e~^2-"^"p da^, (47) 

and 

D_{u' ,u — u') = I W{up^ u — Up) T^{up,u — Up, Up — u') dMp.(48) 

J Up 

4.1.2. Conservation of the wide-field measurement equation 

Appendix I A . 31 demonstrates that the wide-field dirty image 
is here again the convolution of the sky brightness / by a 
wide-field dirty beam D"' or, in the Fourier plane, 

iLy (") ^ (^^){u) = £ !''(«', u - u') I{u') du' (49) 

with 

D^(u', u — u') = / W{up, u — Up) ^ {up, u — Up, u') dup, (50) 

J Up 

where E^(Mp, a^, u') 

= // g(up - Up)7(as - a'J E(Up,a',,Wp - u')dMpda',. (51) 



We thus have equations that resemble those containing the 
sampling function alone, except for 1) the replacement of 
the generalized sampling function E by its gridded version 
and 2) the way the variables are linked together both 
in the gridding of E (i.e., Eq. I5ip and in the averaging of 

(i.e., Eq.Eni). 
4.2. Regular resampling 

It is well known that too low a resampling rate in one 
spa ce implies power aliasing in the conjugate space (See 
e.g. lBracewellll2000l: iPress et al.lll992l) . Aliasing must be 
avoided as much as possible because it folds power outside 
the imaged region back into it. Table |3] defines the inter- 
vals of definition of the different functions we are dealing 
with [i.e., visibilities, primary beam, dirty image, and dirty 
beam), as well as the associated sampling rates needed to 
enforce Nyquist sampling. The boundary values of the def- 
inition intervals (|mL„^ and |q;L,„^) are related to the sam- 
pling rates [da and du, respectively) through 



l"lmax-^"= l"lmax-^W = 



1 



^samp 



(52) 



where risamp is an integer characterizing the sampling. 
Nyquist sampling implies risamp = 2. However, slight over- 
sampling (e.g. risamp = 3) is often recommended because 
the measures suffer from errors and the deconvolution is 
a nonlinear process. In this section, we examine the prop- 
erties of the different functions to define their associated 
sampling rates. 

4.2.1. The Cis sampling rate of the visibility function 

When Fourier transforming the measurement Eq. [T] along 
the as axis, we derive the Ekers & Rots equation p^ . 
This equation implies that V{up,Us) is bounded inside the 
[—rfprim, +(iprim] Spatial frequency interval along the Ug 
axis. As a result, the visibility function needs to be regularly 
resampled at a rate of only 0.5/dpiim to s atisfy the Nyquis t 
theorem. This was first pointed out by ICornwelll 



This sampling rate is equal to 0prim/2 or ~ 0twhm/2.4. The 
"usual, wrong" habit of sampling at ^fwhm /2 is indeed un- 
dersampling with aliasing as a consequence. iMangum et al.l 
(2007) discuss the consequences of undersampling in-depth 
in the framework of single-dish imaging. 

4.2.2. The Up sampling rate of the visibility function 

Now, the Fourier transform of the measurement equation [1] 
along the Up axis gives 



V{ap, Us) = B{ap - as)I{ap), 



where V{ap,as) V(up,as). 



(53) 



(54) 



We use the tilde sign under V to denote the inverse Fourier 
transform of V along its first dimension. A well-known 
Fourier transform property implies that B has infinite sup- 
port because B is bounded. The resamphng rate along the 
Up axis therefore depends on the properties of the product 
of B{ap — as) times I{ap) as a function of Up. While no 
unique answer exists, three facts help us to find the right 
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Table 3. Interval ranges of definition and associated sam- 
pling rates for the used functions. 



Functions Intervals 




Samplings 


Visibilities 


Up 


^ c/max 


dup 


— 2 c/alias /^samp 




Oip 


< ^ali 


as/2 


dap 


— ^syn /^samp 




Us 


< dpr 


m 


dus 


— 2 c/iniagc/''T'samp 




as 


< Qira 


agc/2 


das 


— ^prim / ?^samp 


Primary beam 


u's 


< rfpr 


m 


du's 


— 2 c/aiias /?7-samp 




a^s 


< ^ali 


as/2 


fin,' 

oas 


— f^prim/ '''samp 




It 

Us 


< dpr 


m 


du':^ 


— 2 C^alias/^samp 




< 


1 < Sah 


as/2 


da's 


— ^prim/^samp 


Dirty image 


u\ 






du - 


- 2 c/image/^samp 




a\ 


< ^ima 


gc/2 


da - 


■~ ^syn / ^samp 


Dirty beam 


u 


^ C^max 


du' 


— 2 cZiniagc/^samp 




a' 


^ ^imagc/2 


da' 


— ^syn/?T'samp 




II 

u 


< dpr 


m 


du" 


— 2 c/imagc / ''T'samp 




a" 


1 < 


agc/2 


da" 


— ^prim/^samp 



Table 4. Minimum sizes of the dirty beam images to get 
an image fidelity or a dynamic range greater than a given 
value. 



Minimum fidelity 


^alias 


/ ^fwhm 


a 




or dynamic range 


(/b = 0)'' 


(/b = 


0.0625) 


(/b = 0.1) 


10^ 


2.2 




2.2 


2.2 


10^ 


3.5 




3.7 


6.6 


10* 


8.4 




13.4 


13.7 


10'' 


19.8 




> 20.0 


> 20.0 



Notes. The image sizes are expressed in units of the primary 
beam full width at half maxium. The computation is done 
for 3 different ratios of the secondary-to-primary diameters (i.e. 
/b, the antenna blockage factors). The values are derived from 
the modeling of the antenna power patterns shown in Fig. [3] 



sampling rate: 1) B falls off relatively quickly. 2) The result 
depends on the spatial distribution of the sky brightness 
and in particular on the dynamic range in brightness needed 
to accurately image it. 3) The measure of V{ap,as) has a 
limited accuracy owing to thermal noise, phase noise, and 
other possible systematics (e.g. pointing errors). For sim- 
plicity, we quantify the measurement accuracy by a single 
number, namely the maximum instru mental fideli t y mea - 
sured in the image plane as defined in iPetv et al] ()2001|) . 
There are two cases: 

1. The maximum instrumental fidelity hmits the d y namic 
range in brightness. For instance, iPetv et al] ()20Q1[ ) 
showed that the fidelity of interferometric imaging at 
(sub)-millimeter wavelengths will be limited to a few 
hundred. In this case, V{ap,as) aliasing can be toler- 
ated when the amplitude of B is less than a fraction of 
the inverse of the maximum instrumental fidelity. 

2. The maximum instrumental fidelity is much greater 
than the image fidelity, as can be the case at centimeter 
wavelengths. In this case, V{ap^as) aliasing can only be 
tolerated when the amplitude of B is less than a fraction 
of the inverse of the dynamic range of the image. 

The criterion derived in each case gives a typical image 
size (^aiias), which can be converted into the desired Up 
sampling rate. To be more quantitative. Fig. [3] models the 




5 10 15 20 

0/(0.5 0t^hm) 

Fig. 3. Simple models of the antenna power patterns as 
a function of the sky angle in units of half the primary 
beam FWHM (6 

fwhm)- In the 3 cases shown, the illumi- 
nation is Gaussian with an edge taper of 12.5 dB but 3 
different ratios of the secondary-to-primary diameters {i.e. 
fh, the ani:enna blockage factors) are considered (see e.g. 
lGoldsmithlll998l chapter 6). The middle and bottom pan- 
els respectively model ALMA and PdBI antennas. The red 
lines define the minimum angular sizes for which the an- 
tenna power pattern is less than a given fraction. 



normalized antenna power patterns of an antenna illumi- 
nated by a Gaussian beam of 12.5 dB edge taper and with 
a given blockage factor (ratio of the secondary-to-primary 
diameters). The top panel presents an ideal case with- 
out secondary miror, while the middle and bottom pan- 
els present simple models of the ALMA and PdBI anten- 
nas. The largest angular sizes at which the power patterns 
are less than a given value, 'Pqi is a first-order estimate of 
^alias/2 to get a fidelity or dynamic range higher than l/Vo- 
Table H] gives the values of ^aiias/^'fwhm as a function of the 
searched fidelity or dynamic range. This condition is suffi- 
cient but not necessary. Indeed, the aliasing properties also 
depend on the brightness distribution of the source. 

4.2.3. The u sampling rate of /diity(u) 

We have no garantee that the sky outside the targeted field 
of view is devoid of signal, so the only way to ensure a 
given dynamic range inside the targeted field of view is to 
choose the image size large enough so that the aliasing of 
potential outside sources is negligible. This means that the 



9 



J. Pety and N. Rodriguez-Fernandez: Revisiting the theory of interferometric wide-field synthesis 



dirty image size must be equal to the field-of-view size plus 
the tolerable aliasing size 

6*801(1 + 6'alias- (55) 

The conjugate uv distance and associated uv sampling then 
are 

7 Afield T ^ ^image /rn\ 

dimagc = -; J- — and du = . (56) 

1 I afield ri,„ — „ 



4.2.4. The u' and u" sampling rates of D.{u' ,u") 

The u" axis must thus be sampled at the same rate as 
the second dimension of the definition space of S, i.e., as 
Us. Moreover, u' has in this equation a behavior (u' = 
Up + m") similar to u {— Up + Us). It must thus have the 
same sampling behavior as u. This sampling rate {du' — 
dimage/nsamp) IS Quite high. Some deconvolution methods 
(see below) allow us to relax this sampling rate. 

4.3. Absence of gridding "correction" 

Imaging of single-field observations goes through the fol- 
lowing steps: 1) convolution by a gridding kernel, 2) regular 
resampling, 3) fast Fourier transform, and 4) gridding "cor- 
rection". The so-called gridding "correction" is a division 
of the dirty beam and dirty image by the Fourier transform 
of the gridding kernel used in the initial convolution. This 
step is mandatory when imaging single-field observations 
to keep the image-plane measur ement equation as a sim - 
ple convolution equation (see e.g. lSramek fc Schwablll989h . 
When imaging wide-field observations, as proposed here, 
the Fourier transform along the as dimension, followed by 
the shift-and-average operation, freeze the convolution ker- 
nel into the dirty beam of the wide-field measurement equa- 
tion. This is why the gridding "correction" step is irrelevant 
here. 



5. Dirty beams, weighting, and deconvolution 

In radioastronomy, the dirty beam is the response of the 
interferometer to a point source. In the wide-field synthesis 
framework, the response of the interferometer to a point 
source, D, a priori depends on the source position on the 
sky. D{a',a") can thus be interpreted as a set of dirty 
beams, with each dirty beam referred to by its fixed a" 
sky coordinate. These simple facts raise several questions. 
What are the properties of the convolution kernel? Is it 
possible to modify these properties? How do we deconvolve 
the dirty image? 

5.1. A set of wide-field dirty beams 

With the wide-field synthesis framework proposed here. 
Appendix I A. 61 shows that 
D{a',a") = 

// B{a" — a' — as)^^{a' — ap, a" — as)A{ap, as) dapdas.(57) 

JJ ctpas 

Up 

where A(ap,as) D S{up,as), (58) 



{a' .a") 

and n{a',a") ^^„^ W{u',u"). (59) 

A(ap,as) is the single- field dirty beam, associated with the 
uv sampling at the sky coordinate as- And il,{a' ,a") will be 
called the image plane weighting function, while W{u' ,u") 
is the uv plane weighting function. The set of wide-field 
dirty beams D is then the double convolution of the image 
plane weighting function and the single-field dirty beams, 
apodized by the primary beam at the current sky position 
as. 

While the shape of the single-field dirty beam is directly 
given by the Fourier transform of the sampling function, 
the shape of the wide-field dirty beam depends, directly or 
through Fourier transforms, on the sampling function (5), 
the primary beam shape (i?), and the weighting function 
{W). Moreover, the wide-field dirty beam shape a priori 
varies slowly with the sky position, since it is basically con- 
stant over the primary beamwidth as stated in Sect. 14.21 
It nevertheless varies, implying, for instance, a "slow" vari- 
ation of the synthesized resolution over the whole field of 
view. 

While the single-field and wide-field dirty beam expres- 
sions seem very different, they share the same property of 
expressing the way the interferometer is used to synthe- 
size a telescope of larger diameter in the image plane. In 
other words, the sampling function for single-field imag- 
ing and D. for wide-field imaging express the sensitivity of 
the interferometer to a given spatial frequency. These uv 
space functio ns are called the transfer functions of the in- 
terferometer (|Thompson et al.l ll986. chapter 5). Modifying 
the transfer function has a direct impact on the measured 
quantity. Once the interferometer is designed and the ob- 
servations are done, the only way to change this transfer 
function is data weighting. 

An ideal set of wide-field dirty beams, D{a' ,a"), would 
have the following properties. All the wide-field dirty beams 
should be identical (i.e., independent of the a" sky coordi- 
nate) and equal to a narrow Gaussian (its FWHM giving 
the image resolution). This would give the product of a wide 
Gaussian ofu' by a Dirac function of u" , as the ideal wide- 
field transfer function, Dfu'.u" ). 

5.2. Dirty beam shapes and weighting 

When imaging single-field observations, giving a multiplica- 
tive weight to each visibility sample is an easy way to 
modify the shape of the dirty beam and thus the prop- 
erties of the dirty and deconvolved images. Natural weight- 
ing (which maximizes signal-to- noise ratio) , robust weight- 
ing (which maximizes resolution), and tapering (which en- 
hances brightness sensitivity at the cost of a lower resolu- 
tion) are the most popu lar weighting techniques (see e.g. 
ISramek fc Schwablfigsl . 

In the case of wide-field synthesis, a multiplicative 
weight can also be attributed to each visibility sample be- 
fore any processing. However, the weighting is also at the 
heart of the wide-field synthesis because it is an essential 
part of the shift-and-average operation. No constraint has 
been set on the weighting function up to this point, which 
indicates that the weighting function {W) gives us a de- 
gree of freedom in the imaging process. We look in turn at 
both kinds of weighting. In both cases, an obvious issue is 
the definition of the optimum weighting functions. As in the 
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case of single-field imaging, there is no single answer to this 
question. It depends on the conditions of the observation 
and on the imaging goals. 

5.2.1. Weighting the measured visibilities 

Natural weighting consists of slightly changing the defini- 
tion of the sampling function. It is now set to a normalized 
natural weight where there is a measure and elsewhere. 
The natural weight is usually defined as the inverse of 
the thermal noise variance, computed from the radiometric 
equation, i.e., from the system temperature, the frequency 
resolution, and the integration time. Using this weighting 
scheme before computing the first Fourier transform along 
the as sky dimension makes sense because the observing 
conditions (and thus the noise) vary from visibility to visi- 
bility. 

We propose to generalize this weighting scheme to other 
observing conditions than just the system noise. Indeed, 
critical limitations of interferometric wide-field imaging are 
pointing errors, tracking errors, atmospheric phase noise (in 
the (sub)-millimeter domain), etc. While techniques exist 
for coping with these problem s (e.g., water vapor rad iome- 
ter, direction-dependent gains : [Bhatnagar et al.ll2008[ ). they 
are not perfect. The usual way to deal with the remaining 
problems is to flag the source data based on a priori knowl- 
edge of the problems, e.g., pointing measurement, tracking 
errors, rms phase noise on calibrators, etc. However, flag- 
ging involves the definition of thresholds, while reality is 
never black and white. It can thus be asked whether some 
weighting scheme could be devised to minimize the effect of 
pointing errors, tracking errors or phase noise on the result- 
ing image. We propose to modulate natural weighting based 
on the a priori knowledge of the observing conditions. 

5.2.2. Weighting the synthesized visibilities 

Robust weighting or tapering the measured visibilities do 
not make sense in wide-field synthesis because the dirty 
image is made from the synthesized visibilities after the first 
Fourier transform along the as sky dimension. A weighting 
function W then appears naturally as part of the shift- 
and-average operator. Its optimum value depends on the 
properties of the measured sampling function. Here are a 
few examples. 

Infinite, continuous sampling. This is the ideal case stud- 
ied in Sect. 13.11 Knowing that the Ekers & Rots 
equation (llSp links the quantity we want to estimate, 
i.e., I, to many noisj0 measurements, V{up,Us), via a 
product by B (assumed to be perfectly defined), we 
can invoke a simple least-sq uares argument (see e.g. 
iBevington fc RobinsonI [20031 ) to demonstrate that the 
optimum weighting function is 



W{up, u — Up) 



w{up, u — Up) B{up — u) 
w{up, u — Up) B [up — u) dup 



,(60) 



with w{up,Us) the weight computed from the inverse 
of the noise variance of V{up,Us). Using Eq. [201 it is 
then easy to demonstrate that D{u) = 1, and then 



® The noise is assumed to have a Gaussian probability distri- 
bution function. 



Idi-ctyid) — I (a). The dirty image is a direct estimate of 
the sky brightness; i.e., deconvolution is superfiuous. 

Complete sampling. The signal is Nyquist-sampled, but it 
has a finite support in both the uv and sky planes, im- 
plying a finite synthesized resolution and a finite field 
of view. In contrast to the previous case, this one may 
have practical applications, e.g., observations done with 
ALMA in its compact configuration. Indeed, the large 
number of independent baselines coupled to the design 
of the ALMA compact configuration ensure a complete, 
almost Gaussian, sampling for each snapshot. In this 
case, the best choice may be to choose the weighting 
function so that all the dirty beams are identical to the 
same Gaussian function. In this case, the deconvolution 
would also be superfluous. 

Incomplete sampling. This is the more general case studied 
in Sect. 13.21 The signal not only has a flnite support but 
it also is undersampled (at least in the uv plane). The 
deconvolution is mandatory. The choice of the weighting 
function thus will depend on imaging goals. 
If the user needs the best signal-to-noise ratio, some 
kind of natural weighting will be needed. It is tempting 
to use Eq. [5D| as a natural weighting scheme. However, 
the main condition for derivation of this weighting func- 
tion, i.e., the Ekers & Rots equation (fTS)) . is not valid 
anymore, as the noisy measured quantity {SV) is now 
linked to the quantity we want to estimate (/) by a local 
average (see Eq. I3ip . This is why it was more appropri- 
ate to try to get a Gaussian dirty beam shape in the 
complete sampling case. 

If the signal-to-noise ratio is high enough, the user has 
two choices. Either he/she wants to maximize angular 
resolution power and needs some kind of robust weight- 
ing, or he/she wants to get the more homogeneous dirty 
beam shape over the whole field of view. This require- 
ment cannot always be fully met. The Ekers & Rots 
scheme enables us to recover unmeasured spatial fre- 
quencies only in regions near to measured ones, because 
B has a finite support. 

5.3. Deconvolution 

Writing the image-plane measurement equation in a 
convolution-like way is very interesting because all the de- 
convolution methods developed in the past 30 yea rs are op- 
timiz e d to treat deconvolut i on pr o blems (see e.g. Hogbom 
1975 [Claig Imdt iSchwabI [l985 iNaravan fc Nitvananda 
1986 ^). For instance, it should be possible to deconvolve 
Eq. [23] with just slight modifications to the standard 
CLEAN algorithms. Indeed, Eq. [531 can be interpreted as 
the convolution of the sky brightness by a set of dirty 
beams, so that the only change, once a CLEAN component 
is found, would be the need to find the right dirty beam 
in this set in order to remove the CLEAN component from 
the residual image. 

Following [Clar3 (Il980ll and ISchwabI (11984 . most al- 
gorithms today deconvolve in alternate minor and major 
cycles. During a minor cycle, a solution of the deconvolu- 
tion is sought with a simplified (hence approximate) dirty 
beam. During a major cycle, the current solution is sub- 
tracted either from the original dirty image using the exact 
dirty beam or from the measured visibilities, implying a 
new gridding step. In both cases, the major cycles result in 
greater accuracy. The iteration of minor and major cycles 
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enables one to find an accurate solution with better com- 
puting efhciency. In our case, the approximate dirty beams 
used in the minor cycle could be 1) dirty beams of a much 
smaller size than the image, or 2) a reduced set of dirty 
beams (i.e., guessing that the typical variation sizescale of 
the dirty beams with the sky coordinate is much larger 
than the primary beamwidth), or 3) both simultaneously. 
The model would be subtracted from the original visibili- 
ties before re-imaging at each major cycle. The trade-off is 
between the memory space needed to store a full set of ac- 
curate dirty beams and the time needed to image the data 
at each major step. Some quantitative analysis is needed to 
know how far the dirty beams can be approximated in the 
minor cycle. 

It is worth noting that the accuracy of the deconvolved 
image will be affected by edge effects. Indeed, the dirty 
brightness at the edges of the observed field of view is at- 
tenuated by the primary beam shape. When deconvolving 
these edges, the deconvolved brightness will be less precise, 
because the primary beam has a low amplitude there. This 
only affects the edges, because inside the field of view, every 
sky position should be observed a fraction of the time with 
a primary beam amplitude between 0.5 and 1. This edge 
effect is nevertheless expected to be much less troublesome 
than the inhomogeneous noise level resulting from standard 
mosaicking imaging (see Sect. 17. ip . 

6. Short spacings 

6.1. The missing flux problem 

Radio interferometers are bandpass instruments; i.e., they 
filter out not only the spacings longer than the largest base- 
line length but also the spacings shorter than the shortest 
baseline length, which is typically comparable to the di- 
ameter of the interferometer antennas. In particular, radio 
interferometers do not measure the visibility at the center 
of the uv plane (the so-called "zero spacing"), which is the 
total flux of the source in the measured field of view. 

The lack of short baselines or short spacings has strong 
effects as soon as the size of the source is more than about 
1/3 to 1/2 of the interferometer primary beam. Indeed, 
when the size of the source is small compared to the pri- 
mary beam of the interferometer, the deconvolution algo- 
rithms use, in one way or another, the information of the 
flux at the lowest measured spatial frequencies for extrap- 
olating the total flux of the source. The extreme case is a 
point source at the phase center for which the amplitude 
of all the visibilities is constant and equal to the total flux 
of the source: Extrapolation is then exact. However, the 
larger the size of the source, the worse the extrapolation, 
which then underestimates the total source flux. This is 
the well-known problem of the missing flux that observers 
sometimes note when comparing a source flux measured by 
a mm interferometer with the flux observed with a single- 
dish antenna. 

Wide-fleld synthesis does not recover the full short spac- 
ings. Let us assume that the visibility function is continu- 
ously sampled from d^in to dmax, with d^m ~ l-Sfiprim- 
The length of the averaging hnepatlil, L{u), can be inter- 
preted as the number of measures that contribute to the 
estimation of I{u). Figure |4] shows the variations of L{u) 

^ The notion of averaging linepath has been introduced in 
Sect. 13. Il (see in particular Eg. I16|l . 




u M 

Fig. 4. Length of the averaging linepaths displayed as black 
lines in panel c.2) of Fig.[51 as a function of the spatial scale 
in the flnal, wide-fleld uv plane. In the case of a continuous 
sampling of Up between dmin and dmax, these quantities can 
be interpreted as the number of measures that contribute 
to the estimate of I{u). 



function when starting from a visibility function continu- 
ously deflned in the [c?miii, (^max] interval along the Up di- 
mension. We can expect to recover I{u) only inside the 
[c'min — rfprim, ^max + rfprim] interval. In particular, informa- 
tion on short spacings lower than dmin — dprim (e.g. the 
crucial zero spacing) cannot be recovered when using a 
homegeneous interferometer, and the short spacings in the 
interval [dmin — c^prim, c?min] are recovered with increasing 
accuracy from d^-^ — dprim to dmin- Both effects imply the 
need for complementary instruments to accurately measure 
the missing short-spacings. 

6.2. Usual hardware and software solutions 

To derive the correct result for larger source sizes, it is 
necessary to complement the interferometer data with ad- 
ditional data, which contain the missing short-spacing in- 
formation. The IRAM-30m single-dish telescope is used to 
complement the Plateau de Bure Interferometer. Short- 
spacing information can also be in part recovered with a 
secondary array of smaller antennas and shorter baselines 
(e.g. the CARMA interferometer). In the ALMA project, 
the short-spacing information will be derived by a combina- 
tion of four 12m-single-dish antennas and an interferometer 
of 12 antennas of 7 meters called AC A (Atacama Compact 
Array) . 

From the software point-of-view, two main families of 
algorithms exist in the standard processing of mosaics. 
Either the short-spacing information is combined on the 
deconvolved image (i.e., the interferometer data is im- 
aged and deconvolved separatel y) through a hybr idization 
in the Fourier plane (see e.g. iPetv et all |2001D . or the 
long and short-spacing information is imaged and/or de- 
convolved jointly. In this category, we flnd the pseudo- 
visibility technique, which produces interf erometric-like vis- 
ibilities from single-dish inaps ( see e.g. iPetv et ahl 120011 : 
iRodriguez- Fernandez et ahl I2008L and references therein), 
and the multi-resolution deconvolution algorithms, which 
work on images containing different spatial frequency 
ranges. 

In the next two sections, we show how wide-field synthe- 
sis naturally processes the short-spacing information either 
from single-dish or from heterogeneous arrays. 
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6.3. Processing short spacings from single-dish measurements 

The single-dish measurement equation can be written as 



/sd(a) = Ssd{a) 



/ Bsd{a' -a)I{a')da', 

J a' 



(61) 



where Isd is the measured single-dish intensity, the 
single-dish sampling function, and -Bgd the single-dish an- 
tenna power pattern. As already stated in the introduc- 
tion, the above integral is identical to the ideal measure- 
ment equation of interferometric wide-field imaging taken 
in Up = 0. If we define a single-dish visibility function as 



Vsd{up = 0,a) 



Bsd{a' - a)I{a')da', 



(62) 



we can thus write the measured single-dish intensity as 

Isd{a) = Ssd{a) Vsd{up = 0, a). (63) 

The recognition that the single-dish measurement equa- 
tion is a particular case of the interferometric wide-field 
measurement equation opens the way to treating both the 
single-dish and interferometric data sets through exactly 
the same processing steps. We just have to define a hybrid 
sampling function, 5hyb, as 

Shyh{up ^ 0,a) ^ S{up,a) (64) 

Shyh{up = 0, a) = S'sd(a), (65) 

the Fourier transform of the hybrid primary beam, -Bhyb, 



BhyhiUp ^ 0, u') ^ B{up - u') 



(66) 
(67) 



Bhyh{Up = 0,u') = Bsd{-u'), 
and a hybrid weighting function, IVhyb, as 

WhyhiUp j^O,u' + u" - Up) = Whyb("p,'"' + u" - Up), (68) 

Whyb(% ^0,u' + u") - Wsd{u' + u"). (69) 

All the processing steps described in the previous sections 
(including a potential gridding step of single-dish, on-the- 
fly data) can then be directly applied to the hybrid data 
set. Using the wide- field synthesis formalism, we can easily 
write 



Ihyh{u) = j iZhyb("'> " - ""') ^("') du', 



(70) 

(71) 



with /hyb(w) = -fdirtyl"") -|- M^sd(u) ^sd(w) 

and i2hyb('"'7 

= Riu', u") + Wsdiu' + u")Ssd{u")Bsd{-u'). (72) 
We thus see that /hyb is a linear combination of the infor- 
mation measured by the single-dish (Isd) and by the inter- 
ferometer (/dirty)- There, Wsd(w) plays a particular role for 
two reasons. First, its dependency on the spatial frequency 
(u) enables us to filter out the highest spatial frequencies 
that are measured by the single-dish antenna with low ac- 
curacy. Second, it is well-known that the relative weight of 
the single-dish to interferometric data is a critical parame- 
ter in the pro cessing o f the short s pacings from single-dish 
data (see e.g. iRodrieuez-Fernandez et all[2008i) . This rela- 
tive weight is a free parameter within the restrictions set by 
the noise level {i.e., we want the single-dish data to bring 
information and not just noise to the interferometric data), 
and a criterion must therefore be defined to adjust it to 
an optimal value. We refer the reader to the discussion of 
Sect. [SI which also applies here. 



Table 5. Definition of the symbols used to expose the pro- 
cessing of the short spacings. 



Symbol 


& Definition 


Plane (s) 


/sd 


Measured single-dish intensity 


sky 


Bsd 


Single-dish antenna power pattern 


sky 




Single-dish sampling function 


sky 




Single-dish visibility function 


sky 


W^sd 


Single-dish uii-plane weighting function 


uv 


bi 


Voltage pattern of antenna i 


sky 


Bij 


Power pattern of antenna i and j (= bi b*) 


sky 


V^j 


Visibility between antenna i and j 


uv & sky 


Ihyb 


Hybrid dirty image 


sky 


Bhyh 


Hybrid antenna power pattern 


sky 


Shyb 


Hybrid sampling function 


uv & sky 


Whyh 


Hybrid uv-plane weighting function 


uv & uv 


Dhyh 


Set of hybrid dirty beams 


sky & sky 



Notes. This table uses similar conventions as Table [T] The top 
part defines symbols related to single-dish measurements. The 
middle part defines symbols related to heterogeneous-array mea- 
surements. The bottom part defines hybrid symbols, which re- 
sults from combinations of single-dish and heterogeneous-array 
measurements. 



6.4. Processing short spacings from heterogeneous arrays 

A heterogeneous array is an interferometer composed with 
antennas of different diameters. ALMA and CARMA are 
two such examples. The measurement equation for a het- 
erogeneous array is 



Vij{up,as)= bi{ap-as) b*Aap-as) I{ap) < 



■ dap, (73) 



where bi and bj are the voltage reception patterns of the 
antenna pair that forms the i j baseline and the asterisk de- 
notes the complex conjugate (jThompson et al.lll986i . chap- 
ter 3). The formalism developed in the previous sections 
holds as long as we redefine 

B,j{a) = b,{a)b*ia). (74) 
A simple application of the correlation theorem implies that 



Bij{u) ~ / bi{u + u') bj{u') du' . 



(75) 



The use of the baseline indices ij must be generalized 
throughout the equations because the knowledge of the an- 
tenna type must be attached to each individual data point 
(visibility). As a result, the wide-field synthesis formalism 
can be easily adapted to heterogeneous arrays at the price 
of additional bookkeeping. 



6.5. Two textbool< cases: IRAM-SOm 
ACA 



PdBI and ALMA 



Figure [5] sketches why wide-field synthesis naturally han- 
dles the short spacings in two textbook cases. In the ideal 
case, the Fourier transform along the as dimension pro- 
duces visibilities, which are related to the wide-field spatial 
frequencies of the source brightness weighted by the transfer 
function of the interferometer. In this sense. Fig. [5] displays 
the natural weighting of the synthesized wide-field visibil- 
ities at the position of each measured visibility. Handling 
visibilities from antenna of different sizes just implies that 
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50 100 



uv distance [m] 




50 100 



uv distance [m] 

Fig. 5. Sketches of the natural weighting of the synthesized 
wide-field visibilities. Each measured spatial frequency will 
produce wide-field spatial frequencies apodized by the 
transfer function (B) centered on the measured spatial fre- 
quency. The used transfer function depends on the tele- 
scopes used, explaining why wide-synthesis naturally han- 
dles the short spacing either from a single-dish antenna or 
from a heterogeneous array. The synthesized visibilities in 
the overlapping regions will then be averaged. Two text- 
book examples are illustrated: 1) the combination of data 
from the IRAM-30m single-dish (red transfer function) and 
from the Plateau de Bure Interferometer (black transfer 
functions) at the top; and 2) the combination of data from 
ALMA 12m-antennas used either in single-dish mode (red 
transfer function), in interferometric mode (black transfer 
functions) and of data from the ACA 7m-antennas (blue 
transfer functions) at the bottom. The minimum uv dis- 
tances measured by each interferometer were set from the 
minimum possible distance between antennas (24 m for 
PdBI, 15 m for ALMA and 9 m for ACA). 



the natural weighting function of the synthesized visibilities 
will have a different shape. 

The top panel of Fig. [5] displays how the IRAM-30m 
single-dish is used to complement the Plateau de Bure inter- 
ferometer visibilities. The bottom panel displays how ACA 
is used to produce the short spacing information for ALMA. 
The four 12m- antennas will provide the single-dish informa- 
tion, while the 12 additional 7m-antennas will form with 
ALMA a heterogeneous array. In the first design, ACA and 
ALMA form two independent interferometers; i.e., they are 
not cross-correlated. The singlc-dish antennas, ACA and 
ALMA, thus appear as three different instruments. It is 
thus possible to decompose the hybrid set of wide-field dirty 



beams obtained by processing the 3 sets of data together 
in 3 different sets of dirty beams 

R^y^iu', u") = R,2miu', u'O+^TmK, u")+D^^iu', u"), (76) 

with li^yb = 

Whyh{up, u' + u" - Up) Shyh{up, u") Bhyb{up, u') dup.{77) 
For a multiplying interferometer, 

V|up| < dpi-im, S{up,as)=0 and S{up,Us)^0. (78) 

This implies that Z^sd contributes sd Up — in the sum 

over Up in Eq. [771 Hjm contributes for 9 m < Up ^ 40 m 

and D_i2mi''^' j''^") contributes for 15m < Up < 150m in the 
most compact configuration of ALMA. 

7. Comparison with standard nonlinear mosaicking 

7.1. Mosaicking in a nutshell 

Several excellent descriptions of the mosaicking imag- 
ing and deconv olution algorith ms can be found (see e.g . 
ICornwelll 119881 : ICornwell et all 119931: ISault et al.l [l996bh . 
Here, we summarize the approach implemented in the 
GILDAS/mapping software used to image and deconvolve 
the data from the Plateau de Bure Interferometer. This ap- 
proach is based on original ideas by F. Viallefond in the 
early 90s (|Gueth et al.lil995D . 

The basic ideas of nonlinear mosaicking are 1) imaging 
the different fields of the mosaic independently, 2) linearly 
adding the single- field dirty images into a dirty mosaic, and 
3) jointly deconvolving the dirty mosaic. 

7.1.1. Single-field imaging 

For simplicity, we skip the gridding convolution in the fol- 
lowing equations because the gridding step does not change 
the nature of the equations. Imaging the fields individually 
means that we will work at constant ag . We first define the 
single-field dirty image of the a^-field as 

Op _ 

Isid{oip,as) D Isfd{up,as), (79) 

where the Fourier transform of the single-field dirty image 
is the product of the sampling function S{up,as) and the 
visibility function V{up, as): 

Isfd(up,as) = SV{up,as). (80) 
From the previous equations, it is easily demonstrated that 

/sfd (ckp, as) = / A{ap-ap, as) [B{a'p - as)I{a'p)] da^, (81) 
where the single-field dirty beam is defined as 

Up 

A{ap,as) D S{up,as). (82) 

We can rewrite the previous equation as 

/55(ap) = {A"=*(B"=/)}(ap), (83) 

meaning that the single-field dirty images can be written as 
a local convolution of 5"=/ and A"% the single- field dirty 
beam associated to the currently imaged field. 
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Table 6. Definition of the symbols used to expose the mo- 
saicking formalism. 



Symbol & Definition 


Plane(s) 




Scanned angle 


sky 


Us 


Scanned spatial frequency 


uv 




Phased angle 


sky 


Up 


Phased spatial frequency 


uv 


I 


Sky brightness 


sky 


B 


Primary beam 


sky 


V 


Visibility function 


uv & sky 


S 


Sampling function 


uv & sky 


A 


Set of single-field dirty beams 


sky & sky 


Isfd 


Set of single-field dirty images 


sky & sky 




Set of mosaicked dirty beams 


sky & sky 


^mos 


Mosaicking sky weighting function 


sky & sky 


Wmos 


Mosaicking uv weighting function 


uv & uv 




Mosaicked dirty image 


sky 


mos 


Mosaicked noise image 


sky 


R 


Residual image 


sky 


SNR 


Signal-to-noise ratio image 


sky 



Notes. This table uses similar conventions as Table [T] The top 
part repeats the symbols which have the same meaning in both 
the mosaicking and the wide-field synthesis formalisms, while 
the bottom part defines symbols that either have a different 
meaning in both formalisms or are used only in the mosaicking 
formalism. 



beam and the noise vary across the mosaic field of view. We 
describe here the adaptations made in GILDAS/mapping 
of the simp l est C LEAN deconvohition method, described 
in iHoeboml (|l974l) . Adaptations of more evolved CLEAN 
deconvolution methods are also implemented following the 
same basic rules. 

1. We first initialize the residual and signal-to- noise maps 
from the dirty and noise maps 



Ro(ap) = /mos (ctp) 



and SNRo(ap) 



3(«p) 



(89) 



2. The fc*^ CLEAN component is sought on the SNR^^i 
map instead of the Rk-i map to ensure that noise peaks 
at the edges of the mosaic are not confused with the true 
signal of the same magnitude. 

3. Using that the fc*'' CLEAN component is a point source 
of intensity Ik at position afc, the residual and signal- 
to-noise maps are then upgraded as follows: 

Rkictp) = Rk~i{ap) 



^mos (*^p 1 



-Q!fe, as) B(ak-as) da^, (90) 



7.1.2. Mosaicking the dirty images 

In GiLDAS/MAPPiNcfl, the single-field dirty images are 
formed on the same grid (in particular the same pixel size 
and the same image size covering about twice the mosaic 
field of view). These single-field dirty images are then lin- 
early averaged as 



IraosioLp) = I fltnosiap,as) Isfdiap,as)das, (84) 

w{as) B{ap - as) 



where rimos (ctp, Os) 



(85) 



/^^ w(as) B'^{ap - as) dus 
and w{as) is the sky plane weighting function, i.e., 
u;(as) = ^ (5(q!s - a^)^. (86) 

i * 

In the previous equation, the ai are the positions of each 
sky-plane me asurement, and aj is the rms noise associ- 
ated with Jj^. lCornwell et aTl (|1993[ ) demonstrates that the 
noise in the mosaic image naturally varies across the field 
as 



.^nios i^p) 



(87) 



^ J^^ w{as) B'^{ap - as) das 
In particular, it rises sharply at the edges of the mosaic 



7.1.3. Joint deconvolution 

Standard algorithms of single-field deconvolution must be 
adapted to the mosaicking case because both the dirty 



See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS software. 



and SNRfc(Q!p) 



Rkjctp) 

Nmosiap) 



(91) 



Here 7(~ 0.2) is the usual loop gain that ensures con- 
vergence of the CLEAN algorithms. 
4. Steps 2 and 3 are iterated as long as the stopping crite- 
rion is not met. 



7.1.4. Wide-field measurement equation 

To help the comparison between mosaicking and wide-field 
synthesis, we now go one step further than is usually done 
in the description of mosaicking; i.e., we write the image- 
plane measurement equation as a wide-field measurement 
equation of the same kind as Eq. [531 Substituting Eq. |5l] 
into Eq. [M] and reordering the terms after inverting the 
order of the sum over as and ap, one obtains 

Imosidp) = Dmos{ap - a'p,ap) I{a'p)da'p, (92) 

with Dmos{a',a") 

B{a" -a' - as) n^osia" , as) A(a', a,) da,. (93) 



Taking the inverse Fourier transforms of Umos, we get the 

mosaicking transfer function 

D^osiu',u-u')^ 



Wrr,os{u-Up, Us-u') S{Up, Up-Us) B{up-u') dUpdUs, (94) 



{a' .a") 

with anos(a',a") ^ Wraos{u',u"). 



(95) 
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7.2. Comparison 

While both mosaicking and wide-field synthesis produce 
image-plane measurement equations of the same kind (see 
Eqs.B^land lMI) . the comparison of the dirty beams fEas.[571 
and [93]) and of the transfer functions fEas. l35]and [94|) im- 
mediately shows the different dependencies on the primary 
beams (B), the single- field dirty beams (A), the image- 
plane weighting functions (fl), and their respective Fourier 
transforms (_B, S and W). This means that mosaicking 
is not mathematically equivalent to wide-field synthesis, 
though both methods recover the sky brightness. These dif- 
ferences come directly from the difi'erences in the process- 
ing. If we momentarily forget the gridding steps, mosaicking 
starts with a Fourier transform along the Up dimension of 
the visibility function, and most of the processing thus hap- 
pens in the sky plane, while wide-field synthesis starts with 
a Fourier transform along the as dimension, and most of 
the processing thus happens in the uv plane. 

Moreover, both methods are irreducible to each other. 
Wide-field synthesis gives a more complex dirty beam for- 
mulation in the image plane, which could give the impres- 
sion that it is a generalization of mosaicking. Indeed, the 
wide-field image-plane weighting function can be chosen as 
the product of a Dirac function of a' times a function uj of 
a", 

I.e., n{a',a")^5{a')uj{a"). (96) 

This implies a wide-field uw-plane weighting function in- 
dependent of u'\ i.e., W{u\u") = lJ{u"). This choice is a 
clear limitation because it enables us to infiuence the trans- 
fer function only locally (around each measured Up spatial 
frequency), while weighting is generally intended to globally 
infiuence the transfer function (see Sect. [5]). Eitherway, in 
this case, the wide-field dirty beam can easily be simplified 
to 

D{a\ a") = I B{a" -a'- a^) w(a" - a^) A(a', a^) das.(97) 

J as 

While this simplified formulation of the wide-field dirty 
beam is closer to the mosaicking formulation, they still dif- 
fer in a major way: ui{a"-as) is a shift-invariant function 
contrary to rinios(Q!",Q;s). This is the shift-dependent prop- 
erty of Onios(a",Q^s), which implies the additional complex- 
ity (integral over Ug in addition to the integral over Up) of 
the mosaicking (Eq. [M|) over the wide-field transfer func- 
tion (Eq.[35l). 

One main difference between the two processing meth- 
ods is that standard mosaicking prescribes a precise weight- 
ing function, while we argue that the wide-field weighting 
function should be defined according to the context (see 
Sect. [5]). Another important difference is the treatment of 
the short spacings, which are naturally processed in the 
wide-field synthesis methods, but which needs a very spe- 
cific treatment in mosaicking (see Sect. [6] and references 
therein). Finally, while mosaicking implies a gridding only 
of Up dimension of the measured visibilities, wide-field syn- 
thesis naturally requires a gridding of both the Up and as di- 
mensions. As the Nyquist sampling along the as dimension 
is only 0.5/dprim, the gridding of the sky plane can result 
in a large reduction of the data storage space and cpu pro- 
cessing cost when processing on-the-fiy and/or multi-beam 
observations. 



8. Summary 

Interferometric wide-field imaging implies scanning the sky 
in one way or another (e.g. stop-and-go mosaicking, on-the- 
fly scanning, sampling of the focal plane by multi-beams). 
This produces sampled visibilities SV , which depends both 
on the Mf-plane and sky coo rdinates (e.g., Mh and as). 

Based on a basic idea bv lEkers" fc Rots! (Il979l) . we pro- 
posed a new way to image the interferometric wide-field 
sampled visibilities: SV{up,as). After gridding the mea- 
sured visibilities both in the uv and sky planes, the gridded 
visibilities SV^ are Fourier-transformed along the as sky 

dimension, yielding synthesized visibilities SV sampled on 
a uv grid whose cell size is related to the total field of view; 
i.e., it is much finer than the diameter of the interferometer 
antennas. We thus proposed calling this processing scheme 
"wide-field synthesis" . 

The Fourier transform is performed for each constant 
Up value. As many independent estimates of the uv plane 
are produced as independent values of Up measured. A shift- 
and-average operator is then used to build a final, wide-field 
uv plane, which translates into a wide-field dirty image after 
inverse Fourier transform, i.e., 

^dirtyl"") = / W{Up,U- Up) SV'^ {Up,U~ Up)<lUp, (98) 

J Up 

where is a normalized weighting function. Using these 
tools, we demonstrated that: 

1. The dirty image (7^;^.^^) is a convolution of the sky 
brightness distribution (/) with a set of wide-field dirty 
beams {D^) varying with the sky coordinate a, i.e., 

4rty(«)-/ D^{a-a',a)I{a')da'. (99) 

J a' 

Compared to single-field imaging, the dependency on 
the primary beam is transferred from a product of the 
sky brightness distribution into the definition of the set 
of wide-field dirty beams. 

2. The set of gridded dirty beams (-D^) can be computed 
from the ungridded sampling function {S), the transfer 
function {B, the inverse Fourier transform of the pri- 
mary beam), and the gridding convolution kernel (See 
Eq. EH [501 and EH). 

3. The dependence of the wide-field dirty beams on the 
sky position is slowly-varying, with their shape varying 
on an angular scale typically larger than or equal to the 
primary beamwidth. 

Adaptations of the existing deconvolution algorithms 
should be straightforward. 

A comparison with standard nonlinear mosaicking 
shows that it is not mathematically equivalent to the wide- 
field synthesis proposed here, though both methods do re- 
cover the sky brightness. The main advantages of wide-field 
synthesis over standard nonlinear mosaicking are 

1. Weighting is at the heart of the wide-field synthesis be- 
cause it is an essential part of the shift-and-average op- 
eration. Indeed, not only can a multiplicative weight be 
attributed to each visibility sample before any process- 
ing, but the Mw-plane weighting function {W , see Ea.l98| 
is also a degree of freedom, which should be set accord- 
ing to the conditions of the observation and the imaging 
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goals, e.g. highest signal-to-noise ratio, highest resolu- 
tion, or most uniform resolution over the field of view. 
The W weighting function thus enables us to modify 
the wide-field response of the instrument. On the other 
hand, mosaicking requires a precise weighting function 
in the image plane, which freezes the wide-field response 
of the interferometer. 

2. Wide- field synthesis naturally processes the short spac- 
ings from both single-dish antennas and heterogeneous 
arrays along with the long spacings. Both of them can 
then be jointly deconvolved. 

3. The gridding of the sky plane dimension of the mea- 
sured visibilities, required by the wide-field synthesis, 
may potentially save large amounts of hard-disk space 
and cpu processing power relative to mosaicking when 
handling data sets acquired with the on-the-fly observ- 
ing mode. Wide-field synthesis could thus be particu- 
larly well suited to process on-the-fly observations. 

The wide-field synthesis algorithm is com patible with 
the ity w-unfaceting technique devised by ISault et al.l 
(|1996al) to deal with the celestial projection effect, known as 
non-coplanar baselines (see appendix |B|. Finally, on-the-fly 
observations imply an elongation of the primary beam along 
the scanning direction. These effects can be decreased by 
an increase in the primary beam sampling rate. However, it 
may limit the dynamic range of the image brightness if the 
primary beam sampling rate is too coarse (see appendix [C|). 
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Appendix A: Demonstrations 

A.l. Ekers & Rots scheme 

Fourier-transforming the visibility function along the as di- 
mension at constant Up, we derive with simple replacements 



K,(as)e'^2^"="= da. 



(A.l) 
(A.2) 



B{ap - as) I{ap) e-*2^("p"i>+"»«») da.dap. (A.3) 



We then use the following change of variables /? = ap — a, 
and dl3 — —das, to get 



the theory of interferometric wide-field synthesis 
A.3. G ridding 

The gridding kernel can be defined as the product of two 
functions, each one operating in its own dimension. We use 
this to study separately the effect of gridding in the uv and 
sky planes. We then use the intermediate results to get the 
effect of gridding simultaneously in both planes. 



A. 4. Gridding in the uv plane 

We define the sampled visibility function gridded in the uv 
plane as 



VuJUs) 



(A.4) 



B{(3) I{ap) e-''2- ["p("p + - /^^-l dapd/3 (A.5) 



/(ap)e-*2'^"''("''+"=Map 



S(/3) e-*27r/3(-«,) 

( 

= B{-Us)l{Up + Us). 

A.2. Incomplete sampling 

We here demonstrate that Eg. and Eq.lMlare equivalent. 
To do this, we take the direct Fourier transform of /dirty (a) 



(A.6) 
(A.7) 



I dirty (u) 

D{a - a', a) /(a')e"'^™" dada', 



(A.8) 
(A.9) 



and we replace I {a') by its formulation as a function of its 
Fourier transform 



(A.IO) 



/(u')dw'.(A.ll) 



Using the following change of variables a" = a — a' , a' — 
a — a" and da" = —da' , the innermost integral can be 
written as 



We thus derive /dirty (w) 

Joia - a', Q,)e-*2'^(«"-"'«') dada' 



SVa{up,as) = {g*SV"^}{up) 
g{up-Up)SV°''{u'p)dUp, 



(A.17) 
(A.18) 



Using that the gridding is here applied on the Up dimension, 
while the Fourier transform is applied on the a, dimension, 
it is easy to show that the gridding and Fourier-transform 
operations commute: 

SVlius) (A.19) 

g{up - as)V{u'p, a,)e-^2."="= da3du;(A.20) 



g{up - Up) SVu'^{us) du'p 



(A.21) 



Defining the Fourier transform of the uv gridded dirty im- 
age, we derive 



7^ 

dirty 



{u) = (^SV'yu) (A.22) 
W{up, u — Up) g{up — Up) SVu'^ {u — Up) dMpdUp(A.23) 



Using Eq. |3T] to replace SVu'^{u — Up), we can write the 
Fourier transform of the uv gridded dirty image as 



IdMyiu) = I lf{u',U-u')I{u')Au', 



with 



(A.24) 



{u\u — u') — / W{up,u — Up)^^{up,u — Up,u') dMp(A.25) 



and {up, Us,u') 



D{a - a', a) e-*2vr(a«-a'n') ^^j^/ 



(A.12) 



/ jD{a",a)e 
j ^(u',a)e-'2'^"("-"'Ma 



g{up - Up) S{Up,Us -u -\- Up) B{up - u) du^. (A. 26) 
27ra(u-«')da (A.13) Using ^(Up, Ms - u' + Up) 

~ - . , ^^^^^^ 



D(u\u-u'). 



(A.14) 
(A.15) 



^„.(a3)e-'2."="^ da, 
-9, ^9, 



-i2'K{u —u )cxs 



In the last two steps, we have simply recognized two differ- 
ent steps of Fourier transforms of D. Finally, 



/dirty (") = / I2.{u' ,U - u) I{u')Au' . 



(A.16) 



and S {up,as,u') I) S {up,Us,u'), (A. 28) 

we derive Y,^ {up,as,u') = 

g{up - u'p) as) B{u'p - u') e-'^-K-O-^dw;, (A.29) 
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or A. 5. Gridding in both planes 

f Starting from the definition of SV^ (Eq. E]), we Fourier- 

E (up, as, u') = / g{up - u'^) S(Up, a^, u'^ - u') Au'^. (A. 30) transform it along the sky dimension at constant Up. Using 
"p that the gridding along the Up dimension can be factored 



out of the Fourier transform, we derive 
pling function S. -Qv^ [ n(n,^. _ „M 



Thus, if is the uv gridded version of the generalized sam 



SVl^ius)^ g{up-u'p)SV:,Jus)du'p. (A.42) 



A. 4.1. Gridding in the sky plane 

We defiE 
plane as 



Using Eq. IA.33[ we now replace SVl, (ug) in the previous 
We define the sampled visibility function gridded in the sky equation to get 



SV^Up,as)^{j*SVu^}{as) (A.31) SV^^M ^ j{us) j^^^giup - u'^) SV^,{us) du'^, (A.43) 

7(«s-a'J5K,K)da;, (A.32) 



or Svl^{us)=-f{us)SVl{us). (A.44) 

Applying the convolution theorem on the Fourier transform y^-q^^ this relation, it is easy to deduce that 
along the as dimension, we derive 

^ tF{up,Us,u') ^^{us)^^{up,Us,u'). (A.45) 

SVljus)^^{us)SVu,{us). (A.33) 

Using the convolution theorem when taking the inverse 

Defining the Fourier transform of the sky-gridded dirty im- Fourier transform of along the Ug dimension and re- 



T ' 

dirty V 



age, we derive placing S^(up,a'^,u') with Eq. [QOl we finally derive 

/ ^\ (ur,, as,u') 

Xu)^{sv'^){u) {AM) 

Wiup, u - Up) j{u - Up) SV^^ (u - Up) dup. (A.35) " ^^"^ " < " "'^ Hd«'.-(A.46) 



Using Eq. [JT] to replace SVuj,{u - Up), we can write the A. 6. Wide-field vs single-field dirty beams 

Fourier transform of the sky-gridded dirty image as rm ^ ^- fftti • i i nrr i ii\ 7=:/ / //^ tt • • • 

J to The notation [591 yields W(u . u") = il(u . M l. Usmg this m 

» Eq. [35] gives 

Tlirty(")= / D_\u',u-u')liu')du' (A.36) D{u',u") 

with u-u') " / ^ ~ "P-* '^'^"P' ^^^P ~ ^^P' 

/ — 7 , Taking the inverse Fourier transform along the u" axis of 

Wiup, u~Up)Y, (up,u- Up, Up-u) dup (A.37) Eq. |S3Zland reordering the integral to factor out the term 

independent of u" , we can write 

and ^\up,Us,u'^)=^ius)S{up,Us + u':)Biu':), (A.38) d{u' ,a") ^ [ Biup - u')FTi{up,u' ,a") dup, (A.48) 

or, with the definition of S (i.e.. Eg. [45]). where 

^iup,Us,u'^)=jius)Uup,Us,u';). (A.39) FTi{up,u',a") (A.49) 

_ ^_ =f n{up,u' + u" -Up)Siup,u")e+''^''''"°'" du". (A.50) 

Using Y.''iup,as,u'^) 2^\up,Us,u'^), (A.40) 

We now introduce the following definition 
and the convolution theorem when taking the inverse ^ 

Fourier transform of S^, we derive S{up, u") = I S{up, Ug) e-'27ra.«" ^^^^ (A. 51) 



T7{up,as,u'g) = j{as - a'g)J:{up,a'g,u'^)da'g. (A. 41) 



to derive 
FTi{up,u',a") 



Thus, E is the sky gridded version of the generalized sam- 
pling function E. 



S{up, as) 



n{Up, U' + U" - Mp)e+*27r«"(a"-a,) ^^// 



da.,. 
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Using the following change of variables v = u" + u' — Up, 
u" = V — u' + Up and dv = du" on the innermost integral, 
we get 

FTi(up,w',a") 

S{Up, as)n{Up, a" - as) e+*27r(up-n')(a"-a,) 



B.l. w-axis distortion 

When projection effects are taken into account, the mea- 
surement equation reads 
V{w,Up,as) 



B 



1 -"p 



Substituting this result into Ea. lA.48l and taking the inverse 
Fourier transform along the u' axis, we can write 
D{a',a") = 

Q.{up, a" -as) S{up, as) FT2(up, as,a, a") dupdas, (A. 52) 

JJ UpOis 

where 
FT2(up,as,a',a") 

B{Up - u') e+»2,r«p(a"-a,) g+^2W(a'-a"+a,) ^^f ^ 



In this equation, we continue to work in 1 dimension for 
the sky cosine direction (ap), but we explicitly introduce 
the dependence along the direction perpendicular to the 
sky plane. This dependence appears in two ways, which is 



Using the following change of variables v = Up — u' ^ u' — 
Up — V and dv = du' , we get 



FT2(Mp, as, a', a") = B{a" — a' — ag) e 



-\-i27TUpa 



(A.53) 



Substituting this result into Eq. IA.52I and re-ordering the 
terms, we can write 

L»(a',a")= [ B{a" -a' -as)FT3{as,a',a")das,{AM) 
where 

FT3(as, a', a") = f n{up, a" - as) S{up, as) e+'2-"p"' 

J Up 

A simple application of the convolution theorem gives 
FT3(as, a', a") = / Q{a' - ap, a" - as) A{ap, as) dap, 

J Up 



where A{ap,as) D S(up,as)- 



(A.55) 



Substituting this result into Eq. IA.541 we finally derive the 
desired expression, i.e.. Eg. [571 



Appendix B: From the celestial sphere onto a 
single tangent plane 

Eq. [T] neglects projection effects, known as non-coplanar 
baselines. Any method which deals with interferometric 
wide-field imaging must take this problem into account. 
After a short introduction to the problem, we show how 
wide-field synthesis is compatible with at leas t on e method, 
namely the uvw- unfaceting of ISault et all (fT99 6b'). This 
method tries to build a final wide-field uv plane from differ- 
ent pieces, just as our wide-field synthesis approach does. 
Another promising met hod is the w -projection, based on 
original ideas of Frater fc Dochertvl 111980.') a nd first suc- 
cessfully implemented bv lCornwell et all (|2008t ). We did not 
look yet at its compatibility with wide-field synthesis. 



handled in very different ways. First, the factor y 1 — 
can be absorbed into a generalized sky brightness function 
I{ap) 



I{ap) 



1 - a2 



(B.2) 



After imaging and deconvolution, I{ap) can be easily re- 
stored from the deconvolved I{ap) image. The second de- 
pendence appears as an additional phase, which is written 
as 



P{ap, w) = e-»2.-(Vi^-i). (B.3) 

iThompson et al.l (|1986L chapter 4) shows that this addi- 
tional phase can be neglected only i£| 



Afield 



4 



<C 1 or TT- 



Xd„ 



d^ 

"field 



< 1. 



(B.4) 



The first form of the criterion indicates that the approx- 
imation gets worse at high spatial dynamic range (i.e., 
^'fioid/fi'syn ^ 1) whilc the second form indicates that the 
approximation gets worse at long wavelengths. 

B.2. uvw-unfaceting 

For stop-and-go mosaicking, it is usual to delay-track at 
the center of the primary beam for each pointing/field of 
the mosaic. This phase center is also the natural center of 
projection of each pointing/field. Stop-and-go mosaicking 
thus naturally paves the celestial sphere with as many tan- 
gent planes as there are pointings/fields; i.e., this observ- 
ing scheme is somehow enforcing a uww-faceting scheme. 
In the framework of on-the-f iy observations with ALMA, 
iD'Addario fc EmersonI (|200Cl[ ) indicate that the phase cen- 
ter will be modified between each on-the-fly scan while it 
will stay constant during each on-the-fly scan. This is a 
compromise between loss of coherence and technical possi- 
bilities of the phase-locked loop. Using this hypothesis, the 
maximum sky area covered by the on-the-fly scan must take 
into account the maximum tolerable w-axis distortion. 

The easiest way to deal with such data is to image each 
pointing/fleld around its phase center and then to reproject 
this ima ge onto the mosaic tangent plane as displayed on 
Fig. 5 of lSault et al.l (|1996b[ ). These authors point out that 
this scheme implies a typical u;-axis distortion e less than 



e < (1 — COS^alias) Siu 6'ccntcr ^ TT ^center ^'alias : 



(B.5) 



In contrast to the convention used in this paper, the dmax 
and dfloid unit is meter instead of unit of A in the second form 
of the criterion, in order to explicitely show the dependence on 
the wavelength. 
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where ^center is the angle from the pointing/field center and 
^aiias is the anti- aliasing scale defined in Sect. 14.21 In par- 
ticular, e is at the phase center of each pointing/field. 
In other words, this scheme limits the magnitude of the 
ly-axis distortion to its magnitude on a size equal to the 
anti-aliasing scale (i.e., a few time the primary beamwidth) 
instead of a size equal to the total mosaic field of view. This 
scheme thus solves the projection effect as long as the w- 
axis distortion is negligible at sizes smaller than or equal to 
the anti-aliasing scale. A natural name for this processing 
scheme is uuw-unfaceting because it is the combination of a 
faceting observing mode {i.e., regular change of phase cen- 
ter) and a linear transform of the uv coordinates to derive 
a s ingle sine projection for the whole field of view. 

ISault et al.l (|l996b[ ) also demonstrate that the reprojec- 
tion may be done much more easily and quickly in the uvw 
space before imaging the visibilities because it is then just 
a simple transformation of the uv coordinates, followed by 
a multiplication of the visibilities by a phase term. Finally, 
ISault et al.l (|1996bf ) note that it is the linear character of 
this uv coordinate transform which preserves the measure- 
ment equation [TJ As the change of coordinates happens be- 
fore any other processing, it also conserves all the equations 
derived in the previous sections to implement the wide-field 
synthesis. 



Appendix C: On-the-fly observing mode and 
effective primary beam 

Usual interferometric observing modes (including stop-and- 
go mosaicking) implies that the interferometer antennas ob- 
serve a fixed point of the sky during the integration time. 
Conversely, the on-the-fly observing mode implies that the 
antennas slew on the sky during the integration time. This 
imp hes that the measurement equation [J must be written 
as ([Holdawav fc Fosteij 119941 : iRodriguez-Fernandez et all 

V{up,as) = 

2 rto+St/2 
Jto-St/2 



Table C.l. Definition of the symbols used to explore the 
influence of on-the-fly scanning on the measurement equa- 
tion. 

Symbol & Definition 
St Integration time 
&s Scanned angle averaged during St 
Up Spatial frequency averaged during St 
Sas Angular distance scanned during St 
t^sicw Slew angular velocity of the telescope 

A Primary beam apodizing function 
-Bcfi Effective primary beam resulting 

from OTP scanning: Bcft (a) = {B * A}{a) 
Inearth Angular velocity of a spatial frequency 
due to Earth rotation 



where c?,nax is the maximum baseline length, Woarth is the 
angular velocity of a spatial frequency due to the Earth ro- 
tation (7.27 X 10~^ rads~^), ^aiias and 6'syn are respectively 
the minimum field of view giving a tolerable aliasing and 
the synthesized beam angular values. 



C.2. Effective primary beam 



Assuming that condition IC.3I is ensured, we can write 
Eg. IC. II with the same form as Eq. [T]by the introduction of 
an effective primary beam (-Beff); *-e.. 



V{up,as)= / Seft(ap-d.)/(ap)e-*'""-"''dap, (C.4) 
where i3cff(ap - ds) = — / B{ap ~ as{t)} dt.{C.b) 

"t Jt„-St/2 



rto+St/2 
lto-St/2 

Using the following change of variables 



/3 = a,(i)-d„ dl3^^^^dt or = (C.6) 



we derive 



B{ap - asit)} I{ap) e-^2^«-"-(*) da^ 



dt, (C.l)Seff (ap -as)= B {{ap - d,) - /?} A{I3) d/3 (C.7) 



where 5t is the integration time and Up and d^ are the mean 
spatial frequency and direction cosine, defined as 



^ r-to+St/2 
— / Up{t)dt 
Jta-St/2 



and 



1 

St 



to+St/2 

a,{t)dt.{C.2) 

to-St/2 



with A{I3) 



and Sa^ = 



1 



Wslcw(/3) St 



n 



A 

Sa. 



to+St/2 



to-St/2 



fslcw(i) dt. 



(C.8) 



(C.9) 



In this section, we analyze the consequences of the antenna 
slewing on the accuracy of the wide-field synthesis. 

C.l. Time averaging 

In all interferometric observing modes, it is usual to adjust 
the integration time so that Up{t) can be approximated as 
Up. To do this, it is enough to ensure that Up{t) always 
varies less than the uv distance associated with tolerable 
aliasing (daiias, see Sect. 14. 2p during the integration time 
{St) 



(5i < 



^alias 



^max '^carth 



dt 
1 sec 



6900 



-'alias / "sy 



/ ^syn 



(C.3) 



In these equations, Wsicw(/3) is the slew angular velocity of 
the telescope as a function of the sky position, Sag is the an- 
gular distance covered during St, A is an apodizing function, 
and n(/3) is the usual rectangle function, which reproduces 
the finite character of the time integration. 

C.3. Interpretation 

The form of the measurement equation is conserved when 
averaging the visibility function over a finite integration 
time, as long as the true primary beam is replaced by an 
effective primary beam, which is the convolution of the true 
primary beam by an apodizing function. To go further, it 
is important to return to the two dimensional case. Indeed, 
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Fig. C.l. Assessement of the relative error implied by the 
use of the true primary beam instead of the effective pri- 
mary beam when analyzing interferometric on-the-fly data 
sets. Left: Inverse Fourier transform of interferometer pri- 
mary beam, B (i.e. the autocorrelation of the antenna il- 
lumination). Right: Relative error as a function of sam- 
pling rate of the primary beam. The curves of different col- 
ors show the results at different normalized uv distances 
(u/dprim) from the center of B. 

the convolution must be done along the slewing direction, 
resulting in an effective primary beam elongated in a par- 
ticular direction. 

In principle, the equations derived in Sect. |3] can be ac- 
commodated just by replacing the true primary beam by 
its effective associate. In practice, the probability to take 
into account the effective primary beam is low because its 
shape varies with time. Indeed, it is often assumed that the 
sky is slewed along a straight line at constant angular veloc- 
ity. Even in this simplest case, it is advisable to slew along 
at least two perpendicular directions to average system- 
atic errors, implying two different effective primary beams. 
However, practical reasons may /will lead to complex scan- 
ning patterns: 1) The limitation of the acceleration when 
trying to image a square region leads to spiral or Lissajous 
scanning patterns; 2) The probable absence of derotators 
in future multi-beam receivers (B. Lazareff, private commu- 
nication) implies the need to take into account the Earth 
rotation in the scanning patterns of the off-axis pixels. 



The relative error implied by the use of the true primary 
beam instead of the effective primary beam is then 

B,s{u)-B{u) _1 1 

Bcs{u) A{u) smc{uSas)' 

Fig. IC.ll shows this relative error as a function of the num- 
ber of samples per primary beam FWHM in the image 
plane (i.e., 9iwhm/Sas) for different uv distances (in units of 
dprim)- We see that we derive a 1% accuracy at all u when 
we sample the image plane at a rate of 5 dumps per primary 
beam. However getting a 0.1% accuracy needs quite high 
sampling rates (about 15). This must be compared with the 
accuracy of knowledge of B. 

We note that if a better accuracy is needed than the 
one achievable with the highest sampling rate, it is in the- 
ory possible to replace in the correlator software the rect- 
angle apodizing function by another function which falls 
more smoothly. To avoid the loss of sensitivity inherent to 
the use of such an apodizing function (by throwing away 
data at the edges of the time interval of integration), would 
require, for instance, to half-overlap the integration inter- 
vals. This would imply more book-keeping in the correlator 
software and some noise correlation between the measured 
visibilities. 



C.4. Approximation accuracy 

In the following, we thus ask what is the trade-off accu- 
racy of using the true primary beam instead of the effective 
primary beam. The first point to mention is that using dif- 
ferent scanning patterns somehow helps because the averag- 
ing process then ma kes the bias less systematic. Following 
iHoldawav fc Foster] ((1994) , we quantify the accuracy lost in 
the Fourier plane. Indeed, the Ekers & Rots scheme tries to 
estimate missing sky brightness Fourier components from 
their measurements apodized by the Fourier transform of 
the primary beam. In the Fourier space, the above convolu- 
tion just translates into a product. The Fourier transform 
of the apodizing function thus degrades the sensitivity of 
the measured visibility, V{up,as), to spatial frequencies at 
the edges of the interval [up — dprim , Up + dprim] ■ To guide us 
in our quantification of the accuracy lost, we now explore 
the simplest case of linear scanning at constant velocity, 
where Wsiow(/?) is constant and das = Wsiow (5^. The Fourier 
transform of the apodizing function is then a sine function: 

A{u) ^smc{uSas). (C.IO) 
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